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1. Introduction 

Lanchester-type* differential-equation combat models are an important 
tool for analyzing many important problems of military operations research. 
In such a combat model, a so-called attrition-rate coefficient represents 
the fire effectiveness of a particular weapon-system type against a particular 
target type, i.e. the weapon-system type's effective firepower against such 
a target. Time-dependent attrition-rate coefficients are used to model 
temporal variations in firepower on the battlefield. Thus, we see that 
time-dependent attrition-rate coefficients are important (and, in fact, 
essential [4-6]) for the quantitative analysis of hypothetical combat. 

Militarily realistic computer-based Lanchester-type models of quite 
complex military systems have been developed for almost the entire spectrum 
of combat operations, from combat between battalion-sized units to theater- 
level operations. Nevertheless, a simple combat model may yield a clearer 
understanding of significant interrelationships that are difficult to per- 
ceive in a more complex model, and such insights can subsequently provide 
valuable guidance for more detailed computerized investigations. In this 
report we consider such a simplified variable-coefficient Lanchester-type 
model of combat between two homogeneous forces. 

For this variable-coefficient Lanchester-type model of combat 
between two homogeneous forces, different functional forms for the attrition- 
rate coefficients lead to different mathematical functions being involved 
in representing and computing the force-level trajectories. In a previous 
paper [5] we have discussed the plausibility of the hypothesis that except 


for the special case of a constant ratio of attrition-rate coefficients, 


* 
So-called after pioneering work of F. W. Lanchester [3]. 


ni 


the solutions to such differential equations cannot be represented in term 
of "elementary" functions of analysis. Thus, new transcendental functions 
arise in the study of combat modelled with time-dependent attrition-rate 
coefficients. In particular, we have previously introduced [5-6] so-called 
Lanchester-Clifford-Schlafli (LCS) functions for analyzing combat modelled 
with power attrition-rate coefficients with "no offset" (see Section 3 below). 
In the Appendix to this report is contained a reduced set of tables 
for the LCS functions: it contains tables of five-decimal-place values of 


the hyperbolic-like LCS functions PGs), H 4)» and T Ox) (see Section 4 


1 
below) for 11 fractional values of a (see Section 6 below). A companion 
report [8] contains the most extensive set of tables currently available. 

The main body of this report provides the theoretical and modelling background 
for the use of these tables. In particular, we examine a model of a constant- 
speed attack on a static defensive position and show how associated range- 
dependent kill rates give rise to time-dependent attrition-rate coefficients 
with "no offset.'' Numerical computations are presented to illustrate the 

use of the LCS functions for analyzing such "aimed-fire' combat. As a con- 
sequence of the availability of these tables, one can now study this variable- 


coefficient combat model almost as easily and thoroughly as Lanchester's 


classic constant-coefficient model. 


2. Variable-Coefficient Lanchester-Type Equations of Modern Warfare. 
We consider combat between two homogeneous forces modelled by the 


following variable-coefficient Lanchester-type [3] (see [4,5]) equations 


of modern warfare 


dx 


a -a(t)y with x(0) = XQ > 

. (23) 
Oy (0) = 
ae b(t)x with y(0) Yo» 


where t = 0 denotes the time at which the battle begins, x(t) and y(t) 
denote the numbers of X and Y at time t, and a(t) and b(t) denote 
time-dependent Lanchester attrition-rate coefficients, which represent the 
effectiveness of each side's fire. These coefficients depend on variables 
such as fseee separation, tactical posture of targets, rate of target acqui- 
sition, firing rate, etc. (see [4-7] for further details). Variable attrition- 
rate coefficients are used to model temporal variations in firepower on the 
battlefield. In any analysis of combat, moreover, we should use the above 
equations (2.1) only for x and y > O and, for example, set dx/dt = 0 
when x = 0, since negative force levels have no physical meaning. 
Mathematically, we assume that the attrition-rate coefficients a(t) 
and b(t) are defined, positive, and continuous for to < t < +o with 
t, <0. We also assume that a(t) and b(t) € L(t),T) for any finite 


0 


T > ty. We further take a(t) and b(t) to be given in the form 


a(t) = k a(t) ; and b(t) = kK h(t) : (232) 


where ko and k are positive constants chosen so that a(t)/b(t) = k/k, 
when g(t) = h(t). We introduce the combat-intensity parameter hy and 


the relative-fire-effectiveness parameter Ap defined by 


A, = vk k , and = ks : (2.3) 


R 


From our assumptions about a(t) and b(t), it follows that, for example, 
T 

a(t) f L(t,,T) implies iF a(t)dt = +o. 
0 


The X force level as a function ot time may be represented 


as [5,6] 


x(t) = x, {C, (0) C, (t) _ S,(0)S, (t))} - Yo"rR {c, (0)S, (t) z S,(0)C,(t)}, (2.4) 


where the hyperbolic-like general Lanchester functions (GLF) Cy (t) and Sy (t) 


are linearly-independent solutions to the X force-level equation 


2 
7 1 da 8% _ a(eyb(t)x = 0, (2.5) 


oO 
=~ 
ct 
VY 
i} 
ca) 
. 


Sy (tp) = 0, 
(2.6) 


{1/a(t,)} dc, /dt (to) 


Il 
io) 
w 


{1/a(tp)} ds, /dt (ty) = 1/Vr, . 


Here to denotes the largest finite time at which a(t) or b(t) ceases 
to be defined, positive, or continuous. The Y force level as a function 
_of time is given by a similar expression, with Cy (t) and S(t) being 
analogously defined for the corresponding Y force-level equation. 


It is sometimes convenient to introduce the new independent variable 


t defined by 


E 
T= f Ya(s)b(s) ds . Qs) 
c 
0 


It is readily seen that the transformation 1 = t(t) is well defined and 


invertible. Let us denote 1(0) as To: We observe that Ey = QO implies 


that 1t, > 0. If we denote the "average intensity of combat" as 


0 
Yva(t)b(t) , then 


Seas t 
Ya(t)b(t) t = {(1/t) f Va(s)b(s) ds}t =t- 1 


; (2.8) 
0 0 

The substitution (2.7) transforms (2.5) into 

ane Ai d \ dx 

pe ORre Qn R(t) ac X= 0, (2.9) 

dt 
with initial conditions 

tO Bae ar and {1//R(1,)} dx/dt(t,) = -¥o 


where R(t) = a(t)/b(t). 


3. Combat Modelled with Power Attrition-Rate Coefficients. 

The above equations (2.1) basically apply to "aimed-fire" combat 
when targét-acquisition times do not depend on the numbers of targets available 
(see [5,6] for further details). A large class of tactical situations of 
interest can be modelled with the following general power attrition-rate 


coefficients [5-7] 


a(t) = k(t + oy". and b(t) = k(t +c + A)”, Ga) 


where A and C > 0. We will call A the offset parameter, since it allows 
us to model (with yw and wv > 0) battles between opposing weapon systems 
with different maximum effective ranges (see [5,6]). We will call C the 
starting parameter, since it allows us to model (again, with w and wv > 0) 
battles that begin within the maximum effective ranges of the two opposing 
systems. We observe that for the general power attrition-rate coefficients 
(3.1) we have to = -C, and u and vy must be > -1 in order that a(t) 
and b(t) € L(ty,T). 

The above nomenclature is motivated and possible applications of our 


work are indicated by considering S. Bonder's model of the constant-speed 


attack on a static defensive position (see [4-7] for further details) 


d 
Se ars and rei =B (LH 5 (322) 
where r denotes the range between opposing forces, and a(r) and §8(r) 
denote range-dependent attrition-rate coefficients. Range is related to time 
by 


Ct) = he = Vie (375) 


where Ro denotes the opening range of battle and v > 0O denotes the 
constant attack speed. For example, let us consider the constant-speed attack 


of a homogeneous Y force against the static defensive position of a homo- 


geneous X force. Figure 1 diagrammatically portrays this situation. 
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The basic idea is that force separation, i.e. range between the opposing 
forces, changes over time, and the fire effectiveness of, for example, 
a single Y firer, denoted as a(r), depends on this force separation. 

In many cases of tactical interest, we may model the fire effective- 


ness of, for example, the Y weapon system (as a function of range) with 


a(x) = (3.4) 


where RY denotes the maximum effective range of the Y weapon system and 
u > 0. Here yu is used to model the range dependency of Y's attrition-rate 
coefficient (see Figure 2). We model 8(r) similarly, with corresponding 
quantities Re and v being analogous to Ro and wu _ above. 


If we use (3.3) to eliminate range r from (3.4), we obtain 


dx 
de SES 
(3.5) 
dy _ _ 
cf bCE)x y 


where the time-dependent attrition-rate coefficients a(t) and b(t) are 
given by (3.1). It follows that the offset and starting parameters are given 


-by 
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R, -R R -R 
Le ess , and C = (2a *o), (3.6) 


and that 


u v 
k, = 29() , and ok, = By | : G27) 


a 


We observe that A and C >O if and only if ee R > Ro: By consider- 


ing (3.6) and Figure 3, the reader should have no trouble in understanding 


a 


our terminology for A and C. 


When the offset parameter is equal to zero (i.c. A = 0), then the 


coefficients (3.1) reduce to 


a(t) = k ACHE) and b(t) = k(t +c)” : (3.8) 


We will refer to (3.8) as power attrition-rate coefficients with "no offset." 


As we have seen above in Bonder's constant-speed attack model, these 
coefficients model, for example, combat between weapon systems with the same 
maximum effective range so that there is no "offset" in the "reaching out" 
of the weapon systems against each other in combat (again, see Figure 3). 
For these coefficients (3.8), the transformed X force-level equation (2.9) 


becomes 


Z 
ax + alae dx _ <=... (3.9) 
T dt 


with initial conditions 


10 


[as (Oy?) = 9 Aq usats st AJojowered BuTqzeqs sul “4 

‘as (ra ?a) = y Aq uaat3 st toqzoweied jJasjjo aut -“¢ - (2g yy wnwtu tw > Oy 

(umoyus se) pue 0. se pajousp ST eT33eq Fo asuer Butuado syy “7 *ATaatyvedseaa fy 
pue wy se poqouep aie smaqshks uodeam X pue YX ay} Fo saesuer saTIOaFJa wnutxeu 
OUL ‘T ‘SaLON] ‘HoOeIWe pseds-JueqsUOD BUTTTepow sJUSsTITJJo0O |a}eA-uUOTATIAIIe 


aamod 203 9g dAaeqoweied Butjieqjs oy] pue Yy Aaqeweired Yosjjo oy Jo vuotjeue[tdxy ‘“¢ sean3tTy 


(si9jow) % ebuoy 


OOSE _ 


20 


11 


0) 
$4191} {80D 
@iDY 
“UOT ALY 


9°0 


8°O 


2q-1 


_ ¢ Sere 2a 
x(t) =X > an 5 dt}. Joe oe ee 
0 
Here 
= Vite 
2X 
ye eee (u+v+2)/2 
(Pia) ero , (3.11) 
and 
2X 
_ I (ut+v+2)/2 
to *( geta) A (3622) 


Let us observe that O <q <1 when yw and y > -l. Furthermore, q > 1/2 
if and only if dR/dt <0, i.e. R(t) is a strictly decreasing function of 


time. 


4. Lanchester-Clifford-Schlafli (LCS) Functions. 
Consider the function F, (x) defined by the power series 


w 2k 
F(x) = 1a) J et? (5 
k 


bo fet Met a) r” 


For a #0, -1, -2,... the radius of convergence for F,(x) is infinite 
‘by the ratio test for convergence of power series [2]. Hence, F(z) is 


an entire function of the complex variable z= x + iy, with an essential 


2 


singularity at the point at infinity. Now consider the function HG) 


defined by the infinite series 


©0 (x /2)2 t™) 
H,@) = TG) Lo te racra Fp * het) 
Observing that 
2a 
H (x) = (1/a) (x/2) Fg > (4.3) 


we see that for a> O the infinite series (4.2) is uniformly convergent 
on compact subsets of the complex plane. From (4.3) one can readily deduce 
the recursive relation 


| (ejay 
F Gd) = Foy ® a laa + 1) F gg 6%) : (4.4) 


We will call the functions F Gs) and HG) Lanchester-Clifford-Schlafli 
(LCS) functions (see Note 10 on pp. 66-67 of [5]). Other properties are 
readily deduced and are given in Table I. 

The function F Cx) satisfies the linear second-order ordinary 


differential equation 


Loe eee (4.5) 


with initial conditions 
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Table I. Properties of the LCS Functions F (x) and H (x). 
a a 


bs dF /dx = (x/2)7 7H (x) 


ge dH /dx (x/2)°° Te (x) 


a. For. 1 Gok ay SR 


where @ is not an integer (including zero) 


4. FY Ox=0) = 1 

3% H (x=0) = 0 For J¢.>.0 
6. dF /dx(x=0) = 0 

7 {(x/2)* aH /ax} = 1 

Ss Py) = cosh x 

9. Ay 7) = sirh x 
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dF 
F (0) =l, and a (0) =O , 
while H Go) satisfies 
anu, 2a-L Pes 
2 ( x ) dx as = 0's (4.6) 
dx 
with initial conditions 
x. 1720 dH, | 
H (0) =O, and CG) dx | =], 
) x=0 
| Thus, {Fs Hy J is a fundamental system of solutions to 
| 7 
| 
dF 2a-1l, dF 
a + ( )=-F=0, (4.7) 
x dx 
dx 


with Wronskian W(F )= (fay, It follows that the GLF for the 


H 
» “l-a 


X and Y force-level equations for combat modelled with the attrition-rate 


coefficients (3.8) are given by 


Ay \ 

Cy (t) = Ey) 5, (t) {1 aoey) ; (4.8) 
ry 1-2q 

Cy (t) = Atte) Sy (t) -aies Hey) ; (4.9) 


where p= 1-q. If we define 
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T(x) = Hy GO/EL OO); (4.10) 


then 
S(t) ( hy a H (t(t)) 
1) eo) Vee eo F Gr(t)) ? Met) 
x q 
or 
( dy jee 
Ty (t) = ra a Fatty) ; (4.12) 


where Ty (t) denotes a hyperbolic-like GLF, which corresponds to the 


hyperbolic tangent. Observing that for u, v > -l, lim 44 T(t) = +”, 


we see that TS Ox) is a strictly increasing function of x _ on the interval 


[O, + ©) and 


0<T () < Tew for Ole < Fie, (4.13) 
with 
: T (1-a) 
bine) (4.14) 
Al hy r(a) 


since by the results of Taylor and Comstock [7] the parity-condition parameter 


Q* = Q¥(u, v, C = 0) is given by 


d 2q-1 
Se ees ae ee Tp) 
serie Ty (t) ~ Qk(u, v, 0) (. +vt+ :) T'(q) ° (4.15) 


tr + 


We recall that Taylor and Comstock [7] have introduced the so- 
called parity-condition parameter Q* as the value (or range of such values) 
for the initial condition Q _ to the initial-value problem 


16 


it 
a 
J 


with E(t) 


(4.16) 


It 
ye) 


a Vip b(t)EL with Ey (t 9) 


such that By (t3Q*) and Ey (t3Q*) > 0 for all t 2 to. In other words, 
Q* is the value of Q in (4.16) above such that neither Ey nor E 
ever become zero. In this case, both Ey (t3Q*) and E, (t3Q*) are positive, 
strictly decreasing functions, similar to decreasing exponentials. Thus, 
we may call Q* "the Y equivalent of an X force of unit strength," 
since the forces are "at parity," with neither force being annihilated in 
finite time. Taylor and Comstock have shown that for either a(t) f L(0,+ ~) 
or b(t) £ L(O, + ©), then Q* is unique and given by 

Sy Ct) 


Wo see es 
ve Cy (t) Q* ite 


The significance of the parity-condition parameter Q* is that it allows 


us to predict force annihilation as the following theorem shows. 


THEOREM 1 (Taylor and Comstock [7]): Assume that either 
a(t) f L(O, +2) or b(t) £ L(O, +”). Then the X _ force 
will be annihilated in finite time if and only if 


x C.,(0) - Q*S_(0) 
oe Ne = as (4.18) 


% Q*C, (0) = 8, (0) f 


i 


5. Use of LCS Functions for Analyzing Combat. 
The Lanchester-Clifford-Schlafli (LCS) functions F Ox) and H, GO 


are useful for analyzing "“aimed-fire'" combat (see Section 3 above) modelled 
with the power attrition-rate coefficients with "no offset" (3.8), which 


we rewrite here as 


a(t) = k(t + cyt, and b(t) = k(t + Ne (5.1) 


In other words, the LCS functions arise in solving the differential combat 
model (2.1) with attrition-rate coefficients (5.1). In order that both 
att) and bt) ¢ L(t),T), we must have uw and »v > -l. Military situations 
modelled by these equations have been discussed in Section 3 above, e.g. 
combat between two weapon systems with the same maximum effective range. 
For such combat, the LCS functions may be used to 

(1) compute force-level declines, 

(2) predict force annihilation, 
and (3) predict the time of force annihilation. 

Let us now see how the LCS functions may be used to obtain the 
above information about force-level declines and force-annihilation pre- 
diction. According to (2.4), (4.8), and (4.9) above, the xX force level 


is given by 


x(t) = XQ {F Ct) He ted) - H tg? H(r(t))} 


r 


2q-1 
I 
oe ae a (Fy(tg)H (r(t)) -H (tg) Fo (x(t), 5-2) 
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where q is given by (3.10), p = 1-q, and ;(t) is given by (3.11), which 


we rewrite as 


2) 
Pe ee: yy (HAV+2) /2 
T(t) = ( ae ,) (t + C) ; (5-53) 


The time to annihilate the xX force* is determined by «(t) = 0, and 


it follows that 


| ds \ 
Sopa a" es alae 


¥X 
T(r(t3)) = — ine (5.4) 
xl yo? + Yo 7 eoenred FM 
where from (4.10) 
ee = WE ee ATU) ; (5.5) 


and we recall that p+qe#=I1. It follows that the time to annihilate xX, 


oe is given by 


x 
If we multiply the first equation of (2.1) by y, the second by x, add, 
and integrate, we obtain 


t 


x(t) y(t) = Xg¥q 7 f {a(s) nae + b(s) a te) tae . 
0 


‘ which shows that x(t) and y(t) can have at most one finite zero. Hence, 
af x(tX) = 0, then we know that y(t) > 0 for all t > 0. 
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r q-p 
I 
x. at} | *o'po? * Yo Vip (. ree a: Be oe fo 


= iP 
c t \" 
di 
xg to? 7 a) VR = + 2 ae. 


a q 


Taylor and Comstock [7] have shown that Te is strictly increasing and 


satisfies (see also (4.12) above) 
Se = T'(p)/ tq) 5 . (5.7) 


where p= 1-q. It follows that in order for xX to be annihilated in finite 


time, the right-hand side of (5.4) must be less than [I(p)/I'(q). Let us 


observe that for to == = 0, (5.4) simplifies to 
P-q 
r 
Ty (c(tQ)) = = oe : ee 
Yo Vhp 


Thus, we have proved the following theorem concerning force- 


annihilation prediction. 


THEOREM 2: Consider combat between two homogeneous forces 
modelled by (2.1) with power attrition-rate coefficients 
(5.1). Assume that wu and vy > -l and that the above 
equations hold for all time. Then the X_ force will be 


annihilated in finite time if and only if 
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d q-p 
Pies ees ae a ey) (59) 
| 0 q’ 0 0 rz Lt wt 2 q 0 | , : 

where q = (v+1)/(u + v + 2) and p = 1-q. For ty = 0 
(i.e. C= 0 so that = 0), X will be annihilated in finite 
time if and only if 

X  P(p) ( ‘ ) 

¥ Ta) R\p tv eo 3 en 


6. Tabulation of LCS Functions. 
This report contains a reduced set of tables of the Lanchester-Clifford- 
Schlafli functions. The Appendix contains tables of five-decimal-place values 


of the hyperbolic-like LCS functions F(x), H 9 bX)» and TGs) for various 


4 
values of the argument x, namely x = 0.00 (0.01) 2.00 (0.1) 10.0, and 

a = 1/2, 1/3, 2/3, 1/4, 3/4, 1/5, 2/5, 3/5, 4/5, 3/7, and 4/7. These values 
of the index a correspond to u, v = 0, 1, 2, and 3 in (3.8) and allow one 
to analyze, for example, a basic spectrum of range capabilities for weapon 
systems in the constant-speed-attack model of Section 3. These tables have 
been calculated by the recursive means given in Section 8 of [5]. A more 


extensive tabulation (namely, for a = 1/2, 1/3, 2/3, 1/4, 3/4, 1/5, 2/5, 


3/5, 4/5, 2/7, 3/7, 4/7, 5/7, 4/9, 5/9, 3/11, 5/11, 6/11, 8/11, 5/13, 8/13, 
5/17, 12/17, 5/21, and 16/21 corresponding to u,v = 0, 1/4, 1/2, 1, 1 7 25 3) 


Ze 


is to be found in a companion report [8]. This companion report contains 
the most extensive set of tables of the Lanchester-Clifford-Schlafli 


functions currently available. 


A representative tabulation of the hyperbolic-like LCS functions 
Fix) Hy_,&); and T Gd) is given in, for example, Tables 8A and 8B of 
the Appendix for a = 3/5. The values of the argument x are the same as those 
used for the tabulation of the hyperbolic functions by Abramowitz and Stegun 
[1]. We observe from Table 8B and (4.13) that the limiting value of 
TO) as x ++ (here a = 3/5) is quickly reached, with three-decimal- 
place accuracy already attained for x= 4.5. Moreover, the use of these 
tables (specifically, Tables 8A and 8B of the Appendix) for combat analysis 


is illustrated in the next section. 


7. Numerical Examples 

In this section we examine a couple of numerical examples to show 
some of the insights that may be gained into the dynamics of combat between 
two homogeneous forces from our results (see also [6]). These examples 
illustrate the use of the LCS functions Fix), Hy > and TO) for 
analyzing "aimed-fire'" combat modelled with the power attrition-rate coefficients 
with "no offset" (5.1). As in [4-7], we consider S. Bonder's model (3.2) 
for the constant-speed attack against a static defensive position. We will 
focus on the use of the LCS functions for predicting force annihilation, 
since the computing of force-level trajectories with Lanchester functions 


is adequately handled elsewhere (see [4-5]), 


ze 


Let us accordingly consider the constant-speed attack of a homogeneous 
Y force against the static defensive position of a homogeneous X _ force 
(see Section 3 above tie further modelling details, especially Figure 1). 
For our numerical computations, we assume that the fire effectiveness of the 


Y weapon system varies linearly with range, i.e. 


a(r) = (7.2) 


and that the fire effectiveness of the X weapon system varies quadratically 


with range, i.e. 


MD 
oS 
a 
ray 
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|r 
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B(r) = (7.2) 


with RY = Ra» i.e. both weapon systems have the same maximum effective range. 
In other words, p = 1 in (3.4) and v= 2 for 8(r). We consider a battle 
modelled by the input data given in Table II. In terms of time as the 
independent variable, the attrition-rate coefficients (7.1) and (7.2) become 


via (3.3) 


a(t) = k(t + C) add b(t) = k(t + a (7.3) 
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Table II. Input Data for Numerical Examples 


0.06 X casualties/minute/Y firer 


ao = 
By = 0.6 Y casualties/minute/X firer 
Ro = Re = 2000 meters 

v = 5 miles/hour 
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where Ry= Ra, 


RY - Ry _ Gav z 
Cc = ———__ , k = : and k, = Ata . (7.4) 


WD 


From the input data given in Table II, we compute the parameter values shown 


in Table III, since the transformed X force-level equation is given by (3.9) 


with q = (v+1)/(v+v+ 2), p = l-q, »p = 1, and v = 2. Thus, the xX force 
level may be computed with Pt) and Hy oe with a= q = 3/5. Force- 
annihilation prediction involves the limiting value of T (Tt) = Hy Ct) /F (1) 

a -a a 


as t ++, From Table 8B of the Appendix and Table III, we note the predicted 
agreement between P(1-0)/P (a) and the limiting value of Ty Ox) as. x * 7 = 
[recall (4.13)] for a= q = 3/5. We now consider two cases: (T) Ry = 2000 
meters, and (II) Rg = 1250 meters. 

When Ro = 2000 meters (see Figure 3 of [4]), we have C€ = 0 and 


et O. The maximum time that the battle can last is Soe Ro/v = 14.91 minutes, 
since at this time the attackers reach their final objective, i.e. the 

defender's position (again, see Figure 1). We now consider the qualitative 
behavior of the wp =1, v= 2 force-level trajectory shown in Figure 3 of 


[4]. Theorem 2 tells us that the X force can be annihilated in finite time 


if and only if 


Xp r(p) Ay q-p 
La OO 
Yo Pq) Ke (ees) : (73) 


where q = 3/5 and p = 1-q. Using the numerical values in Table III, we 
compute from (7.3) that the X force can be annihilated in finite time if 
and only if 
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Table III. Parameter Values for Numerical Examples 


kK = 4.0233 
ky = 2.6979 
P= 215,04 
T'(p)/T(q) = 
A=0 


x 


1 


107° xX casualties/(minute)"/y firer 


107° Y casualties/ (minute) ’/X firer 


-48951 
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=— = (0.420 5 (7.4) 


When the X_ force can be annihilated, its annihilation time is given by 


(5.8), which we rewrite here as 


x XQ ae os 
T (t(t_)) = a ; C725) 
q a aa por wok 2 
Yo “Ap 
where 
ie | (utvt2) /2 
7 Ut : 
T(t) = oe ; (7.6) 


Thus, for the numerical values given in Table III, the time of annihilation 


of the X force is given by 


x 
iGo) 454... (7.7) 
q a YQ 
with q = 3/5. We will now illustrate further computations for Xo = 10 and 


Yo = 30. From (7.4) we see that the xX force can be annihilated in finite 


. : X : 
time (but we must verify that t. he A In this case (7.7) becomes 


1.18192 . (7.8) 


T (1(t*)) 
q a 


We must now determine t(t2) such that re) = T (1.18122) by using inter- 


polation methods and Tables 8A and 8B. From Table 8A, we find 
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T = 1, 18172 for = 1,01 
re = 1.17630 for = 1.00 
so that using linear interpolation, we obtain 
X 
T(t.) = 1,009 ; G79) 
whence use of (7.6) yields 
X , 
t= 14.24 minutes, C7210) 


which is less than — = 14.91 minutes so that the defending xX force is 
indeed annihilated before the attacking Y force reaches its final objective. 


Since r(t) = R, - vt, we find that force separation at the instant of 


0 


annihilation of the xX force is 


r ee ae ema): 


Further results may be computed in a similar fashion and are given in 


Table IV. 
When Ry = 1250 meters (see Figure 3 of [5]), we have C = 5.5923 
minutes, tT, = 0.0975, and t = R./v = 9.32 minutes. In this case 
0 max 0 


Theorem 2 tells us that the xX force can be annihilated in finite time 


if and only if 
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Table IV. Annihilation of the xX Force as a Function 


of the Initial Force Ratio for Ro = 2000 meters 
(x,/y,) * (minutes) rs (meters) 

eee i ae ao 
0.333 14.24 89.8 
0.250 11.61 443.2 
0.200 10.19 633.2 
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with q = 3/5 and p =1-q. Using linear interpolation, we obtain from 


Tables 7A and 8A of the Appendix that for the numerical values of Table III 


rt = 2006.2 4 rr = 0.044 , 
C7213) 
F(t.) = 1.0045 Ht.) = 0,223.5 
q 0 pO 
so that (7.12) says that the X force can be annihilated if and only if 
ao 
== < (),382 (7.14) 


When the X force can be annihilated, its annihilation time is given by 


(5.4), which we rewrite here as 


- dy P-q 
a (esl Bete, po 
X ig See See ea eee ar Se ON te on Me ed 
Pg, - , (7.15) 


Pot 
r 
ee ae H(t) 
q 0 tv + 2 q O 


whence for the data of Table III 


3.565u, + Osteo 


X 
Ta‘t(t)) = O156u) + 1.008” ve 
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where Uy = Xo/Yo- Let us also record here that (3.11) yields 


7 (fey tale Cc. (7.17) 


2/ (utv+2) 
a 


We will again illustrate further computations for Xo = 10 and Yo = Bus 

From (7.14) we see that the xX force can be annihilated in finite time (but 
X 

again we must investigate whether or not a cok ne In this case (7.1€) 


becomes 
x 
pone = £33651 4 (7.18) 


whence Table 8A of the Appendix and linear interpolation yield 


ee ae? (7.19) 
so that by (7.17) 
X ‘ 
t= 10.63 minutes . (7.20) 
Since t = R./v = 9.32 minutes < r, we see that the attacking Y force 
max 0 a 


overruns the defender's position before annihilation of the xX force occurs. 


Thus, the battle ends with x, = x(t ¢) > 0 and we? O at t.=t 


f max 


9.32 minutes. Corresponding to te = 9.32 minutes is tT, = 1.1318, and 


then Table 8A of the Appendix yields 


31 


aT 22.138) = 2.909. ee = 1973 4 (a2) 


whence via (2.4), (4.8), (4.9), and (7.13) we obtain 
oS ele.) = ele = 0) =. 25 (7.22) 


Some further numerical results are given in Table V, Again, these parametric 
results should be contrasted with the single’ » =1, v= 2 _ force-level 


trajectory shown in Figure 3 of [5]. 


8. Final Remarks 

In the previous section above, we have seen how the LCS functions 
allow one to conveniently obtain much valuable information about the model 
(2.1) with power attrition-rate coefficients (3.8) without having to explicitly 
compute the entire force-level trajectories. Previously we were limited to 
computing only force-level trajectories (see [4-5]). With the availability 
of these tabulations of LCS functions (see the Appendix of this report and 
[8])}, we can now tell who is going to be annihilated and when this event will 
happen without having to compute the trajectories. Not only did we answer 
questions about the qualitative behavior of the model (e.g. force annihilation) 
for specific values of, for example, initial force levels but also for a 
range of values of the initial force ratio (i.e. parametric analysis of 


model behavior). 


2 


Table V. Annihilation of the X Force as a Function 


of the Initial Force Ratio for Ro = 1250 meters 


(x9/'Yq) tS (minutes) r lereees 
0,333 10.63 bs 
0.250 7.56 235.9 
0.200 6.17 422.8 

te = 9.32 minutes and x, = x(r=0) = 1.35. 

max f 
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The results of this report may be used for other parametric 
analyses, e.g. parametric dependence of battle outcome on attrition-rate 
coefficients. Thus, eis aontenes of this report allow one to develop 
important insights into the dynamics of combat between two homogeneous 
forces with temporal variations in fire effectiveness. With the avail- 
ability of tabulations of the LCS functions, one can now analyze such 
combat modelled by the power attrition-rate coefficients (3.8) with some- 
what the same facility as he can for the constant-coefficient case and thus 
aid in parametric analyses. For further discussions of the significance 


of such results for military operations research, the reader is directed 


to [6] and [7]. 


34 


1s 


REFERENCES 


M. ABRAMOWITZ and I. STEGUN (Editors), Handbook of Mathematical 
Functions, National Bureau of Standards Applied Mathematics Series, 
No. 55, Washington, D.C., 1964. 


K. KNOPP, Infinite Sequences and Series, Dover, New York, 1956. 


F. W. LANCHESTER, "Aircraft in Warfare: The Dawn of the Fourth Arm- 
No. V., The Principle of Concentration," Engineering 98, 422-423 (1914) 
(reprinted on pp. 2138-2148 of The World of Mathematics, J. Newman 
(Editor), Simon and Schuster, New York, 1956). 


J. G. TAYLOR, "Solving Lanchester-Type Equations for ‘Modern Warfare’ 
with Variable Cofficients," Opns. Res. 22, 756-770 (1974). 


J. G. TAYLOR and G. G. BROWN, "Canonical Methods in the Solution of 
Variable-Coefficient Lanchester-Type Equations of Modern Warfare," 


Opns. Res. 24, 44-69 (1976). 


J. G. TAYLOR and G. G. BROWN, "Further Canonical Methods in the 
Solution of Variable-Coefficient Lanchester-Type Equations of Modern 
Warfare,'' NPS55-77-27, Naval Postgraduate School, Monterey, California, 
June 1977. 


J. G. TAYLOR and C. COMSTOCK, "Force-Annihilation Conditions for 
Variable-Coefficient Lanchester-Type Equations of Modern Warfare," 
Naval Res. Log. Quart. 24, 349-371 (1977). 


J. G. TAYLOR and G. G. BROWN, "A Table of Lanchester-Clifford- 
Schldfli Functions," NPS55-77-39, Naval Postgraduate School, Monterey, 
California, October 1977. 


35 


APPENDIX: Tabulation of the LCS Functions F(x), er (x), and T (x) for 
=O Qa 


Oe 1/25 1/35 2/35 2/43/46 5.1/55° 2/55 3/5504) 3563/7 and 4) oy 
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